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Abstract. Floating potential structures that are correlated over a long distance are 
observed with a 2D probe array in the plasma edge of the TJ-II stellarator. We 
introduce a method based on the Singular Value Decomposition to extract the spatio- 
temporal structure of the global, fluctuating, zonal-flow-like floating potential from the 
combined measurements of a 2D probe array and a distant single probe. The amplitude 
of these global structures is seen to modulate not only the high kg spectral power of 
the local turbulence, but also particle transport into the unconfined Scraped-Off Layer, 
as observed by H a monitors around the device. These observations provide the first 
direct evidence of the global modulation of transport by zonal flow-like structures. 
The ability to identify spontaneous and collective rotation events with flux surface 
symmetry opens up the possibility to perform unperturbative studies of the effective 
viscosity in stellarators and tokamaks. 



PACS numbers: 52.25. Os,52.30.-q,52.55.Hc 
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1. Introduction 

Mass flows are an active research topic in magnetic confinement fusion due to their 
importance for plasma stability and confinement. In recent years, both experimental 
and theoretical efforts were made to improve understanding of the momentum transport 
mechanisms that determine the observed plasma rotation profiles (see, e.g., [1]). 
It is generally accepted that fluctuating zonal flows (ZF) in magnetically confined 
plasmas can be driven by turbulence through an inverse cascade process familiar 
from 2D hydrodynamic turbulence j2]. These structures are thought to regulate 
turbulent transport by a) draining energy from small scale turbulence, via the turbulent 
momentum flux (Reynolds and Maxwell stresses), and depositing it into the kg = 
mode (with Eg = —dg(f> = and thus v r = E e /B = ), and b) breaking up turbulent 
eddies once the associated sheared flow pattern is formed. In this view, large-scale flows 
(mean and fluctuating) and small-scale turbulence form a dynamical system whose 
equilibrium state determines the turbulent transport levels [3]. Possibly, ZFs could 
explain one of the most fundamental, long-standing issues in fusion physics: the Low- 
to-High confinement transition (see [4] for a recent experimental finding in this respect). 
Indeed, the practical importance of zonal flows stems from their capability to affect 
global cross-field transport. 

Experimentally, global, fluctuating flow patterns manifest themselves as electric 
potential structures exhibiting long-range correlation (LRC), and have been observed 
in different devices using Langmuir probes and Heavy Ion Beam probes (5j EJ [7J [8]. 
These LRC potential structures oscillate but are constant on a flux surface (which 
amounts to having fcn ~ k± ~ 0), but to qualify as a zonal flow-like structure they 
must also be 'zonal' (i.e. have a finite radial scale k r ^ 0) and be dominated by low 
frequencies (somewhere between the characteristic time-scales of turbulence, ~ 10/is, 
and confinement, ~ 10ms). The turbulent drive of the ZF has mainly been investigated 
[3 [101 [IT] using bi-spectral techniques [I2j . and a flow generation mechanism by 
three-wave coupling was identified, in accordance with rather general considerations 
concerning the structure of drift-wave equations [3]. Although some devices have 
reported a reduction of the broadband spectral power of turbulence associated with 
the presence of large scale flows (see [TBI E], [H] and references therein), the crucial 
practical property of a ZF (i.e., its ability to regulate global transport) has not been 
tested experimentally in fully developed broadband turbulence. 

In this paper, we provide empirical evidence of global transport regulation by 
zonal flow-like potential structures with long distance correlation. For that purpose, we 
present a simple method to isolate the temporal evolution of the fluctuating, collective 
floating potential oscillations from the combined measurements of a 2D Langmuir probe 
array and a toroidally separated single probe pin. This allows us to characterize the time- 
resolved dynamics and to study its effect on edge potential profiles, local turbulence, and 
global particle transport. We find that the instantaneous ZF amplitude modulates global 
particle fluxes dynamically, as reflected by H a radiation monitors around the device. We 
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Figure 1. Pin arrangement in the 2D probe and orientation with respect to the 
magnetic field structure of TJ-II. 



show that the ability to identify spontaneous and collective rotation events with flux 
surface symmetry opens up the possibility to study effective viscosity in stellarators and 
tokamaks, without the need of external momentum input through NBI [15J or external 
biasing [H] [17] . 



2. Experimental Set-up 

Experiments where carried out in the 4-period, flexible, low shear stellarator TJ-II using 
electron cyclotron heated plasmas. The typical magnetic field strength is ~ 1 T and 
the rotational transform j-(a) = 1.65 at the last closed flux surface (LCFS) (typical 
minor radius a is 20 cm), dropping slowly to a central value of j-(0) = 1.55. Electron 
temperatures range from 1 — 0.8 keV at the centre to an edge value of 100 — 50 eV. The 
floating potential measurements presented in this work were obtained with a 2D array of 
5x4 tungsten pins, with radial and poloidal separations of 6 and 3 mm, respectively (see 
figured]). Simultaneous measurements are taken with another pin in a different sector of 
the machine, having a long parallel connection length with the 2D probe (> 100 m).The 
sampling rate of the floating potential signals is 2 MHz. 



3. Spatio-temporal characterization of long range correlated structures by 
SVD 

The results presented here depend crucially on the separation of the ZF contribution 
to the floating potential oscillations -the long range correlated part of the potential 
fluctuations - from the uncorrelated, local fluctuations. This is achieved by means of a 
Singular Value Decomposition (SVD), described below (an alternative method based on 
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the same idea but using the Hilbert-Huang transform has recently been proposed 
early applications of SVD to fluctuations in plasmas can be found in [19]). 
Given a m x n matrix M, its SVD consists of a product of three matrices 

M = UYV^ , (1) 

where U mXm and V nxn are unitary matrices (UW = VV^ = I) and S mxn is a 
real, non-negative, diagonal matrix with elements a a sorted in decreassing order, i.e. 
( <7i > er 2 > . . .). The U and V column vectors, {u a } a=1 . m and {v a } a=1 n are 
orthonormal vector sets and are called topos and chronos respectively [20]. In terms of 
the topos and chronos vectors, equation (Tj[|) can be written as 

M = ^a Q u Q (v a )t . (2) 

a 

This nomenclature is most appropriate when the matrix M consists of the signals 
of spatially separated measurements. This is also the case here, i.e. we construct 
the matrix M as M ik = where the {0/j}j=i...2o are the 20 floating potentials 

measured with the 2D probe and <pf2i is the floating potential of the distant probe. The 
k index ranges over the number of samples in the floating potential signals (6 x 10 4 in 
the dataset used below). 

The cross correlation matrix of the signals for time lag can be expressed in terms 
of the topos vectors: 

n m 

c* = E MhWfjitk) = £ ^<«r • (3) 

k=l a=l 

This can be seen from 

n 



dj = £ M lk M* k ^C = MM^ 



k=l 



Inserting the SVD decomposition (Eq.([T])) one obtains 

C = MM^ = UTy ] (UYVi) ] = UEE T tf , 

which is the matrix form of equation ()3]). The total signal energy has the simple 
expression 

m n 

i=l k=l a 

With our naming convention of the floating potentials, the long range correlation 
between any one probe i = {1 . . . 20} in the 2D array and the remote probe is given by 
Ci 2i- Thus, according to Eq. ()3]), one can write the LRC as a sum of the contributions 
of each of the topos a 

^ c = £^Ki)x = 5>f , (4) 

a a 

we have defined a new quantity, the Partial Eigenmode Correlation z a = cr^(n2 1 )*u a . 
In view of the fact that this quantity gives the correlation of the distant probe with 
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Figure 2. Correlation and frequency analysis of the first two SVD modes of a 30 
ms record of floating potentials. Left column: The z-value of the surface plot shows 
the Partial Eigenmodc Correlation (zf , see text) while the surface coloring represents 
the total long-range correlation (Cf RC = J2 a z f) an d is the same in the two graphs. 
Central column: chronos vectors. Right column: normalized frequency spectrum of 
the chronos vectors (red circles indicate the weighted mean frequency. To ease their 
comparison the spectra of the two chronos are plotted in both graphs, the blue curve 
being the spectrum corresponding to the chrono on the left). 

all the 20 signals in the 2D probe, this vector has 20 entries. The relation (jlj) simply 
expresses that for a particular topo a to contribute significantly to the LRC, it must 
carry a significant fraction of the signal energy (er„) and involve the distant probe signal 

(«&)■ 

The surface plots shown in figure |2] show the Partial Eigenmode Correlation of the 
first two SVD modes of a 30 ms record of floating potential measurements. The average 
plasma density is stable at around (n e ) = 0.57 x 10 19 m -3 , and the innermost pins of 
the 2D probe display a long range correlation of 0.53 (normalized to [—1, 1]). The 20 
first components of z a corresponding to the 2D probe array are rearranged in a square 
matrix, mimicking their actual spatial distribution. The surface colouring corresponds 
to the total LRC of the 2D pins that is given by the sum of the individual Eigenmode 
Partial Correlations according to Eq.(jH). 

This empirical decomposition reveals three characteristics features that are 
consistently reproduced for other similar shots having LRCs. 

• The topo (and the LRC) displays a symmetry in the poloidal direction. While 
this is to be expected for sufficiently low kg modes it has to be stressed that this 
symmetry is not imposed by the SVD technique but appears naturally. 

• The first mode cr 1 u 1 (v 1 )^" accounts for most of the LRC observed in the 2D array 
(compare the colorbar and the ^-values of the topos, i.e. the Partial Eigenmode 
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Figure 3. Temporal evolution of one of the 2D probe floating potentials, the distant 
floating potential, and the first SVD chrono, showing an obvious mutual correlation. 

Correlation, in figEJ). 

• Ther first mode accounts for a significant part of both the cumulative signal energy 
E = Yl™=i °a (namely o\ =1 jE = 58%) and the distant probe's fluctuation energy 

E21 = ELi \Mt k )\ 2 = £™ =1 ^KI 2 (with oLil«rT/3n = 42%). 

• The frequency spectrum of the first chrono is dominated by low frequencies < 10 
kHz. 

From these observations, we conclude that the coherent floating potential evolution given 
by the first SVD mode is a global, low-frequency-dominated, low- fee electric potential 
structure that we will refer to as the Zonal Flow (ZF) component. In matrix form 

<p ZF = a 1 u 1 (vy (5) 

Therefore, the first chrono will be tentatively identified with the evolution of zonal flow 
amplitude, symbolically AzF{tk) — i v l)*- Subtraction of the ZF component from the 
total 2D floating potential evolution provides a local, largely uncorrelated turbulent 
component given by 

01ocal = 4>f ~ 4>ZF (6) 

Both ZF and local components are extracted at the same time resolution as the original 
signals. Also, one should keep in mind that the identification made in equation [5] is not 
absolute in any sense (hence: zonal flow-like), but rather, motivated by the observations 
made above. The conditional averages and correlation analyses presented in the next 
section confirm this interpretation. Figure [3] shows the ZF amplitude Azf{1) (i-e. the 
first SVD chronos) alongside one of the 2D-array floating potentials and the distant 
probe signal. AzF{t) clearly reflects most of the temporal evolution of either potential. 



4. Global modulations of floating potential profile and ZF dynamics 



The modulation of the floating potential profiles 4>f(r] 
electric fields E r = — and shearing rates ojexb = 



by the ZF and their derived 
are shown in figure HI 



1 d 2 4> 



Dynamic 
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Figure 4. Averaged radial profile of floating potentials measured with the 2D probe 
for different ranges of ZF amplitude. The inset shows the probability density function 
of the ZF amplitude. Instantaneous profiles falling within each of the nine bars are 
averaged and correspond to the entries in the legend, where a is the standard deviation 
of the ZF amplitude. The top plot shows the corresponding E r (curves, left axis) and 
shearing rates (circles, right axis). 



Time slices are selected corresponding to specific ranges of Azf values, and the floating 
potentials are time-averaged over those slices. The inset in figure H] shows the probability 
distribution of Azf whose values are given in the normalised units used to represent 
the chronos vectors. The floating potential modulations show how this value translate 
into Volts at the different locations. Thus, we find that positive zonal flow amplitudes 
increase the E x B shearing rate, reaching values over 2 x 10 6 s" 1 at p ~ 0.92, whereas 
negative amplitudes act to flatten the mean potential profile, reducing the shearing rate. 
The overall mean profile coincides with the black curve, corresponding to the zero ZF 
amplitude. The maximum ZF amplitude is located around p = 0.87 near the second 
innermost row. This is confirmed by measurements in which the probe was located 
further inside the plasma column. From the profile modification, a rough estimate of 
the ZF radial scale is ~ 2 cm ^ 20 x (ion Larmor radius). The oscillating radial 
electric field computed from the floating potential gradient varies from 0.5 to 8 kV/m 
at the location of the peak E r . The E x B poloidal velocities obtained from this E r 
estimate range from 0.6 to 10 km/s. 
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Figure 5. Top: Conditionally averaged floating potentials (the conditional trigger is 
Azf/p > 1) showing a similar evolution of the 2D-array and distant probes (average 
values have been subtracted showing differences w.r.t. the mean potential profile - cf. 
figurc[4]). Bottom: Dynamical parameter of a ZF burst extracted from the conditionally 
averaged Azf pulse shape and the model given by Equation[7] Azf units are arbitrary. 
The peak amplitude (6 x 10~ 3 ) correspond to about 20V increase in the inner floating 
potentials. 



It should be noted that these estimates of E T and ujexb do not take into account 
the electron temperature contribution to the plasma potential. Besides, all but the 
innermost pins in the 2D probe are shadowed by the insulating Boron Nitride probe 
body, which may introduce further (although smaller) corrections to the mean plasma 
potential caused by the density wake of the parallel ion flow past the probe [21]. However 
their observed modulations or relative changes are less sensitive to these corrections and 
are considered a good proxy of plasma potential modifications. 

These oscillations are collective for they are observed in both probes. This is 
illustrated by the floating potential evolution shown in figure El conditionally averaged 
over 152 events with A Z f/°~ > 1- The response of the distant probe has the same 
waveform as that observed with the 2D array. The average waveform of the zonal flow 
burst can be used to estimate characteristic parameters of the ZF dynamics. We model 
its time evolution with the simple ID differential equation 

jA ZF (t) = R(t) - i*A ZF (t) (7) 

where R(t) is the ZF drive (which may depend on Azf) and fi is the damping rate 
due to constant viscous forces. An estimate of fi is obtained from an exponential fit 
to the tail of the ZF decay in the [30, 120] /xs interval (figure |5]), assuming the flow is 
freely decaying (R « 0). The resulting decay time is r = (20 ± 2)/is, similar to the 
decay time measured in earlier work [T7] using external biasing switch-off experiments. 
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A ZF 

Figure 6. Expected values of power in the high kg part of the spectrum as a function 
of the ZF amplitude. The absolute spectral power in these wavelengths increases for 
decreasing radius. To simplify comparison, we have normalized the spectral power 
at each radius to the corresponding time average power, and expressed the result as 
percent difference w.r.t. this value. 

Note that this damping rate is a measure of effective viscosity corresponding to the flux 
surface. Using this value, we also obtain the time evolution of the forcing term R(t) by 
integration of Eq.((7]). Figure [5] shows that the waveform of R if almost square. Rise and 
decay times are both ~ 3yus for a total boost duration of 24/is. In physical units, the 
incremental force (acceleration) and asymptotic velocity are respectively -R(O) = 2 x 10 8 
m/s 2 , and R{0)/fi = 4 km/s. 

Equation ([7]) can be interpreted to model the evolution of the flux surface averaged 
poloidal velocity (see, e.g., [3]), where R is the turbulent drive (associated with the 
Reynolds stress) and the —fi term accounts for the neoclassical viscosity or magnetic 
pumping effect. To obtain the dynamical model parameters we had to rely on the 
hypothesis that the forcing term is small in the region of the fit, as suggested by the 
exponential waveform typical of freely decaying flows (see [17], [H]), and assume that 
the viscous coefficient is constant in time. Future work will investigate both the relation 
between the ZF drive and the locally measured Reynolds stress, and the relation between 
the observed damping time and the effective neoclassical flow decay rate. 

5. Modulation of local high-A^ turbulent spectrum and effects on global 
particle fluxes 

The 4 poloidally separated pins of the 2D probe allow us to obtain the instantaneous 
spectral power in the kg = 2h/\q « 0,5 and 10 cm -1 wave- vectors. In figure we 
plot the expected value of the spectral power of the uncorrelated part of the floating 
potential in the two non-zero k e wave numbers, P k > = -P(5 cm -i) + -P(iocm- 1 ); as a function 
of the ZF amplitude for several radial locations. (That is, E(P k > \Azf) where the spatial 
Fourier transform is applied to the corresponding entries of the local floating potential 
</>iocai only, as defined in Eqs. ([5]) and ([6])). It is found that positive ZF amplitudes are 
correlated with a ~ 15% reduction in the small scale fluctuations, while a comparable 
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Figure 7. Cross correlation curves between the zonal flow amplitude and different 
H a monitors around the TJ-II stellarator (inset, blue line). TJ-II has four equivalent 
quadrants, labelled A to D, and each quadrant is subdivided into 8 sectors. The 
H a signals are labelled according to their quadrant and sector. The cross correlation 
between the uncorrelated part of a floating potential (refer to equation [6]) and the H a 
signals is low (inset, grey line). 



increase is observed for negative ZF amplitudes. This is consistent with the profile 
modifications expected in either situation (figHJ), as positive ZF amplitudes increase the 
shearing rates, while negative ZF amplitudes decrease them. 

The effect of these collective radial electric field modulations in the global particle 
transport becomes apparent in figure [7] where the ZF amplitude is plotted together with 
different monitors of H a radiation distributed around the device. H a is mainly emitted 
by neutrals recycled from the walls when outward particle fluxes reach the limiters and 
other material vessel components. Clearly, specific ZF relaxation events are followed by 
a burst of H a radiation in all monitors. This prominent visual impression is confirmed 
by the cross correlation of the zonal flow amplitude and an H a trace shown in the inset. 
The observed 100 /xs time lag in the correlation curve implies an effective radial velocity 
of 100 m/s, assuming that the plasma-wall distance is ~1 cm. This velocity matches the 
typical radial velocity obtained from the turbulent E x B convective transport. On the 
same axis, we also plot the cross correlation of an H a signal with the uncorrelated part 
of the floating potential measured at one of the innermost pins, showing a low positive 
correlation level. Overall, these results suggest that the extraction method presented 
above is successful in identifying ZFs that modify radial transport. 

The results shown here focus on the detailed analysis of one shot, however similar 
behaviour (collective profile modulations anticorrelated with H a radiation) is observed 
for other shots with a similar density (around ~ 0.6 x 10 19 m~ 3 ) displaying long range 



Dynamic transport regulation by zonal flow-like structures in the TJ-II stellarator 11 



correlations. In previous work [22], it was shown that LRCs in TJ-II emerge in the 
vicinity of this value (where the mean electric field changes from a positive to a negative 
root) and at high densities (~ 2x 10 19 m~ 3 ) close to the L-H transition. The investigation 
of the spatio-temporal characteristics of the correlated structures at high density is left 
to future work. 



6. Conclusions 



We have introduced a general method to extract the spatio-temporal evolution of global 
flow structures from the time traces of several 'nearby' floating potential measurements 
and a 'distant' one. This collective mode, dominated by low frequencies (and zonal 
flow-like) is seen to modulate the high kg component of the turbulent spectrum and, 
for the first time, has been shown also to modulate global outward particle fluxes, as 
monitored by various H Q detectors around the TJ-II stellarator. The ability to identify 
spontaneous and collective rotation events with flux surface symmetry opens up the 
possibility to perform unperturbative studies of the effective viscosity in stellarators 
and tokamaks. 
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